An electrostatically denned serial triple quantum dot charged with few electrons 
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A serial triple quantum dot (TQD) electrostatically denned in a GaAs/AlGaAs heterostructure is 

■ characterized by using a nearby quantum point contact as charge detector. Ground state stability 
diagrams demonstrate control in the regime of few electrons charging the TQD. An electrostatic 
model is developed to determine the ground state charge configurations of the TQD. Numerical 
calculations are compared with experimental results. In addition, the tunneling conductance through 
all three quantum dots in series is studied. Quantum cellular automata processes are identified, 
which are where charge reconfiguration between two dots occurs in response to the addition of an 

. electron in the third dot. 

■ PACS numbers: 73.21.La, 73.23.Hk, 73.63.Kv, 81.07.Ta 
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I. INTRODUCTION 

Extensive experimental work has recently been aimed 
towards electrostatically defining and controlling semi- 
conductor quantum dots (QDs) and double quantum dots 
(DQD) iii2iMi§i£i!2i2iI2iii The complete control of the QD 
charge, down to the limit of only one trapped conduction 
band electron, has been demonstrated by monitoring the 
single electron tunneling current through the device 2 -, by 
counting the charge on the QD electron by electron by 
means of a nearby quantum point contact (QPC)^ or 
by combining both methods^. Such efforts are predomi- 
nantly motivated by the desire to control and understand 
the physics of quantum systems, and provide impetus for 
proposals for using the spin! 2 -, charge states^ 3 - or encoded 
subspace a 14 ' 15 of localized electrons as qubits, the ele- 
mentary registers of the hypothetical quantum computer. 
Recent experiments have demonstrated the realization 
and coherent control of charge^ and spin qubits^^ in 
DQDs. 

Extending (double) QD circuits towards a few electron 
triple quantum dot (TQD) is a natural step towards scal- 
able multi-qubit systems. In addition, the spin states in 
three tunnel-coupled QDs can be used to encode a qubit 
in the logical states of a decoherence-free subspace. In 
this way the coherence time of the qubit is expected to 
increase and gate operations to be simplified at the cost 
of a higher number of required QDsi 14 ' 16 d 7 

Quantum information processing relies on a coupling 
between qubits allowing coherent exchange of quantum 
information 1 " 2 -. In QD-based implementations of quan- 
tum computing, where qubit coupling is local, introduc- 
ing coherent qubit transport is important in the design of 
a scalable fault-tolerant architecture.— Coherent trans- 
fer by adiabatic passage (CTAP) has been proposed as a 
way to efficiently move electrons along chains of tunnel 



coupled QDs and entangle quantum mechanical states of 
distant qubitsj 19 i 20 A TQD is the smallest system that 
in principle allows the implementation of CTAP. 

In addition to applications in quantum information 
processing, the interest in TQDs is triggered by a rich 
spectrum of phenomena going beyond the physics of 
DQDs. These include combined charging and reconfig- 
uration events that can be identified as quantum cellu- 
lar automata (QCA) processes^ applications as current 
rectificrS ) 22 i 23 creation of spin-entangled electronS ) 24 i 25 
and new aspects of the Kondo 26 ' 27 i 28 and Fano 2 ^ effects. 

Several efforts have been undertaken to produce later- 
ally defined TQDs. In an early attempt, large TQDs in a 
serial configuration were studied via transport measure- 
ments as a function of the coupling between the QDs. 30 : 31 
More recently, current rectification effects were observed 
in devices consisting of three tunnel coupled QDs charged 
with many electrons. 22 ' 23 ^ 2 Charge stability diagrams at 
low electron numbers were first investigated in a geom- 
etry in which one of two coupled QDs is split further, 
thus realizing a TQD in a ring-like device consisting of 
three tunnel coupled QDsi^ 3 . The mapping of charge sta- 
bility diagrams revealed a QCA effect near points of res- 
onant transport. Magneto-conductance experiments fur- 
ther unveiled Aharonov-Bohm like oscillations^ The re- 
alization of three laterally coupled vertical quantum dots 
is under investigation^ 

In this article we report on the realization of a TQD 
in a serial configuration charged with few electrons. The 
gate layout was specifically designed to define three small 
QDs tunnel coupled in series. We characterize the TQD 
by means of stability diagrams. Integrated charge de- 
tection is performed using a nearby QPC. In addition, 
electron tunneling transport through the three QDs in se- 
ries is investigated. Observed features, which are specific 
for a TQD, including triple points, quadruple points, and 
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QCA effects, arc discussed in detail. Wc derive a classical 
electrostatic model that allows us to predict charge sta- 
bility diagrams of a TQD by minimizing its free energy. 
The model is easily scalable to larger systems containing 
more QDs. A detailed comparison between this model 
and our data is presented in a regime close to points of 
sequential resonant transport through the TQD. 



case of DQDs£i£ The described approach using later- 
ally defined surface gates is in principle scalable to much 
larger systems containing many QDs. 



III. ELECTROSTATIC MODEL OF A TRIPLE 
QUANTUM DOT STABILITY DIAGRAM 



II. THE TRIPLE QUANTUM DOT LAYOUT 

Our sample is fabricated from an AlGaAs/GaAs het- 
erostructurc with a two-dimensional electron system 
(2DES) embedded 120 nm below the surface. At T = 
4.2 K the 2DES features an electron sheet density of 
n s ~ 1.8 x 10 15 m~ 2 and a mobility of \x ~ 75m 2 /Vs. 
Experiments are performed in a dilution refrigerator at 
an electron temperature of T2DES — 100 mK, as deter- 
mined by the width of Coulomb blockade conductance 
peaks<2£ 

Electron beam lithography is used to produce Ti/Au 
gates on the surface of the heterostructure as shown in 
Fig. [TJ The TQD and up to three QPCs are defined by 




Figure 1: (Color online) SEM micrograph of the sample struc- 
ture. Gate electrodes (bright tone) are used to electrostati- 
cally define a TQD and three QPCs. The approximate po- 
sition of the three QDs (A, B, and C) are depicted by black 
circles. Large (small) arrows mark possible tunneling current 
paths through QPCs (the TQD). Ohmic contacts are labeled 
with roman numbers. The gates marked with a, /3, and 7 are 
used as plunger gates of the three QDs A, B, and C. Gates 
marked as dl, d2, and d3 serve to define QPCs as charge 
sensors. 

applying appropriate negative voltages to the gates to lo- 
cally deplete the 2DES beneath. The gate layout extends 
a single QD geometry that allows transport spectroscopy 
at low electron numbers^ Our sample allows the defini- 
tion of up to three QDs ( A, B, and C as indicated by 
black circles in Fig. [1} tunnel coupled in a serial configu- 
ration. Transport measurements can be performed even 
in the regime of only few electrons charging the QDs. 
Three independent QPCs (marked with arrows in Fig.[T]) 
can be used to determine the charge configuration of the 
TQD in the same way as has been demonstrated for the 



A charge stability diagram affords a quick and intuitive 
mechanism to understand many of the properties of a 
quantum electronic system. As the surface gate voltages 
are varied, the system tries to minimize its free energy by 
exchanging electrons with the leads and by redistribut- 
ing the charges between its constituents. In the case of 
a DQD the stable charge configurations form a charac- 
teristic honeycomb diagram as a function of the voltages 
applied to two plunger gates^i The TQD introduces fur- 
ther complexity and richness of phenomena. The obvious 
choice for a full description of all possible charge con- 
figurations of a TQD would be a three-dimensional sta- 
bility diagram as a function of three plunger gate volt- 
ages. Here, we investigate two-dimensional slices of such 
a three-dimensional stability diagram. 

Standard electrostatic models describing a DQD 3 -^ can 
be extended towards a TQD i 32 i 37 We introduce a scal- 
able matrix approach describing electrostatic Coulomb- 
interaction by capacitance matrices. Quantum mechani- 
cal tunneling between QDs is not taken explicitly into 
account for the classical model. Implicitly, tunneling 
of electrons allows transitions between charge configu- 
rations. Fig. [2] sketches an equivalent circuit diagram 
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Figure 2: Equivalent circuit diagram for three tunnel coupled 
QDs A, B, and C in serial configuration. QDs A and C are, 
in addition, tunnel coupled to leads II and III, respectively. 
Tunnel barriers are modeled as resistors and capacitors in 
parallel and electrostatic coupling to three plunger gates a, 
13, and 7 as capacitors. 

for a serial TQD. It models tunnel barriers of the TQD 
as ohmic resistors and capacitors in parallel and electro- 
static coupling to three plunger gates as capacitors. Dur- 
ing a typical measurement all other gate voltages are kept 
constant. They are not included in the equivalent circuit 
diagram in Fig. [5] for simplicity. The relevant circuit con- 
sists of charge nodes (QDs A, B, and C), voltage nodes 
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(plunger gates V a ,Vp : and Vy), and capacitors separating 
nodes. The electrostatic potential of the 2DES including 
source and drain leads is assumed to be at ground level 
(i.e. Vji = Viu = 0). This is a good approximation for 
a typical measurement in the linear response regime. 

In a structure with N nodes with electrostatic poten- 
tials Vj , we can express the total charge Qj of each node 
(including voltage nodes) as 



N 

fe=i 



JV 



(Vi-V k ), 



k=l 



where and Cjk are the continuous polarization charge 
and capacitance between nodes j and k. This expres- 
sion is more conveniently written in matrix notation as 
Q = £V, where Q and V are vectors with the ele- 
ments Qj and Vj, respectively. The matrix £ contains 
the capacitances Cjk between nodes. The diagonal ma- 
trix elements Cjj = Cf are the self capacitances of nodes 
j, defined as the sum of the capacitances between the 
node and all other objects in the universe. Explicitly, 
our model only considers variable voltages applied to the 
plunger gates. All other gates have constant electric po- 
tentials but contribute implicitly via the self capacitances 
Cf. 

We can separate the matrix equation Q = £V as 



(1) 



where charge nodes (QDs) arc labeled with subscript 
D, and voltage nodes (gates) with V. The vectors 

Qd = [Qa,Qb,Qc] , Qv = [Qa,Q/3,Q~f] , V D = 

[Va,Vb,Vc] , and Vy = [V a , Vp, V 7 ] T contain the total 
charges and voltages on the three QDs and three gates, 
respectively. The matrix C is split into 
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-Cab 
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-Cai 
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EC 



-Cac 
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C^ 



containing only capacitances between QDs and the self 
capacitances of the QDs, 



~CAa — Cap — Ca~{ 

-Cbo —Cb(3 — CB7 
-Cca —Cq/3 — Cc 7 



and its transpose containing all capacitances between 
gates and QDs, and Cyv containing all capacitances be- 
tween the three plunger gates. 

In our experiments the electrostatic potentials on all 
gates are independent of the capacitances between the 
gates, because they are imposed by externally applied 
voltages. Hence, the matrix £yv only influences the zero 
point of energy. For simplicity we assume £vv = with- 
out loss of generality. 

Our intention is, to find the ground state stability di- 
agram of a TQD by numerically minimizing its free en- 
ergy F — U — W . Here, U is the electrostatic energy of 



a given configuration created by achieving the work W. 
It is useful to introduce the total effective charge of the 
QDs defined as the sum of Qd and the electrostatic in- 
fluence charge -CdvVv- The relevant part of eq. Q]theii 
reads 



Q D ff = £ D dV d 



Qd - CdvVv 
eN D - CdvVv, 



where we allow only discrete values of the charges of the 
QDs expressed by Qd = eNo with the electronic charge 
e. The vector Nd = [Na, Nb, Nc] t contains the number 
of electrons per QD and defines the charge configuration 
(N A , N b ,N c ) of the TQD. The free energy reads 
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u-w 
1 

2 

\ (Cdd^Qd^ 



V d 
V v 

T 

Qd 



2^0f (E A Qf + E A bQb S + EacQI?) 
— Q% s (EabQ\ S + EbQ% S + EbcQ'c) 
^Q° c ff (ExcQf + E BC Qf + E c Qf) , (2) 



where 
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and X, Y, Z stands for the cyclic permutations of A, B, C. 
In accordance with references and [32| we define the pre- 
factors Ea, Eb, and Ec in eq. as charging energies 
of the individual QDs and Eab, Ebc, and Eac as the 
electrostatic interdot coupling energies between two QDs. 

The charging energies, electrostatic coupling energies, 
and capacitances in cq. ([2]) can be obtained from mea- 
surements, i. e. charge stability diagrams, and the con- 
ductance of the TQD in the non-linear regime. Equa- 
tion @ only takes voltages explicitly into account, that 
are applied to the plunger gates (V a , Va, and Vy). It does 
not consider the detailed geometry of the TQD-device. 
Therefore, the free energy of a given configuration of the 
TQD is not completely determined by eq. Selection 
of a suitable zero point of the charge distribution scales 
a modeled stability diagram to fit measured data. This 
zero point might be defined as the charge on each QD at 
grounded plunger gates. The described model is suited 
rather for a qualitative than a quantitative analysis. 

Figure [3] shows a model stability diagram of a serial 
TQD calculated with eq. as explained above. Here, 
the y and x axes correspond to the plunger gate voltages 
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Figure 3: (Color online) Numerically calculated ground state 
stability diagram of the TQD device shown in Fig.[T]for charg- 
ing energies, electrostatic interdot coupling energies, and ca- 
pacitances between QDs and plunger gates, similar to experi- 
mentally derived values. The color scale of the lines is chosen 
to simulate a possible measurement of the transconductance 
of the left QPC in Fig. [T] as a function of the plunger gate 
voltages Vol and V y , where V~, is modulated. The background 
color denotes zero transconductance. Stable charge configu- 
rations are labeled by triples of numbers (Na, Nb, Nc)- 

V a and V 1 of the two outer QDs (compare Fig. [T]). All 
other gate voltages are kept constant. Lines mark borders 
of stable charge configurations (Na, Nb, Nc)- 

The variable brightness of the lines in Fig. [3J simulates 
an experimental situation, where the left QPC in Fig. [1] 
would be used as charge detector. The brightness reflects 
the electrostatic coupling strength between the QDs and 
the left QPC. In a corresponding experiment the change 
of the current through the QPC in response to an ampli- 
tude modulation of the plunger gate voltage V 1 will be 
measured. Adding the charge of one electron to the TQD 
while increasing V 1 decreases the current through the de- 
tector QPC. This results in a charging line of negative 
transconductance d/Qpc/dV^. 

The spacing between charging lines belonging to a QD 
is approximately proportional to the charging energy Ex 
of that QDj^ The slope of a charging line is always neg- 
ative and determined by the ratio of the respective cou- 
plings between a QD and the two plunger gates a and 
7. Accordingly, the stability diagram of a TQD con- 
tains charging lines with three different main slopes. The 
nearly vertical charging lines indicate charging events of 
QD C, which couples strongly to plunger gate 7 but 
weakly to a (compare the gate layout in Fig. [1]) . Lines 
corresponding to charging events of QD B (in the center) 
have a slope of dV a /dVy = — 1, since we assumed equal 
capacitances between QD B and the two plunger gates a 
and 7. The predetermined symmetry properties result in 
charging lines with reciprocal slopes for QD A compared 
to C. 

All crossings of two charging lines are avoided, be- 
cause of the electrostatic interdot couplings. The result 
are pairs of two triple points, each with three degener- 
ate charge configurations. If well separated from charg- 
ing lines of the third QD, the distance between the two 



triple points of a pair is proportional to the correspond- 
ing electrostatic interdot coupling energy Exy- Since 
our model neglects quantum mechanical tunnel couplings 
these avoided crossings are of purely classical nature. The 
triple points of a pair are connected via a charge reconfig- 
uration Une£& Along these lines with positive slopes an 
increase of Vj always causes a charge transfer between 
the three QDs with the center of charge moving away 
from the detector QPC, hence, resulting in a positive 
transconductance. The total charge of the TQD stays 
constant at the charge reconfiguration lines. 

Between lines of extremal transconductance the 
ground state charge configuration is stable and, hence, 
the transconductance is zero. The electrostatic interdot 
couplings can lead to the zig-zag course of charging and 
charge reconfiguration lines as clearly observable in the 
range of approximately six electron charges on the TQD. 
More complex behavior, including QCA-processes, is ex- 
pected, where charging lines of all three slopes are close 
by, as will be discussed in section V. 

The model stability diagram in Fig. [3J shows the situa- 
tion expected for a TQD in the few electron regime. The 
lack of charging lines in the lower left corner of the figure 
indicates that here, the TQD is uncharged. Along the 
horizontal (vertical) axis QD C (A) is charged electron 
by electron. The plunger gate voltage Vp is chosen such, 
that QD B can only be charged, if another QD is already 
occupied by at least one electron. However, increasing Vp 
would shift the charging lines with slope dV a /dV 7 = — 1 
of QD B towards the lower left corner of the figure. The 
QDs A and C are separated by QD B and hence, have 
a relatively small mutual interdot coupling. This results 
in pairs of triple points being close to each other and 
charging lines that almost intersect (compare e. g. the 
transition between configurations (0,0,1) «-> (1,0,0) in 
Fig. [3]). In comparison, the electrostatic interdot cou- 
plings between neighboring QDs is much stronger result- 
ing in a larger distance between triple points (e. g. see 
charge reconfiguration line between (0,2,4) (1,1,4) 
in Fig. [3J. Note, that for the discussed model calcula- 
tion we chose the coupling between QDs B and C to be 
smaller than that between QDs A and B. 

Resonant tunneling transport of electrons through 
the TQD is only possible at quadruple points, where 
four charge configurations are degenerate. However, 
as quadruple points arc distinct points in a three- 
dimensional space, two-dimensional stability diagrams of 
a TQD containing quadruple points are rare. Since two 
charging lines can meet (but never cross) in one point 
of a stability diagram, a quadruple point of a TQD al- 
ways represents a meeting point of two charging lines and 
two charge reconfiguration lines. A charge stability dia- 
gram in the direct vicinity of quadruple points contains 
up to eight triple points at four avoided crossings. A de- 
tailed discussion of this complex situation and compar- 
isons with measured stability diagrams follows in section 
V. 

In the case of a high degree of symmetry, i. e. equal in- 
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tcrdot couplings Eab = Ebc = ^CA, very rare hextuple 
points with six degenerate charge configurations are the- 
oretically possible. Hextuple points involve the meeting 
of two charging lines and the crossing of two reconfigu- 
ration lines in one point. However, in our serial TQD 
geometry, where two electrostatic intcrdot coupling en- 
ergies arc larger than the third one, we would not expect 
to see such hextuple points. 

Two important limits restrict the validity of the elec- 
trostatic model. The geometry of the electronic probabil- 
ity distribution inside a realistic TQD lacks perfect sym- 
metry. It rather is a complicated function of applied gate 
voltages and the local disorder potential. This causes 
a non-linear gate voltage dependence of the capacitance 
matrix elements. It can result in a change of the distance 
between parallel charging lines or even in a continuous 
change of slopes of charging lines. In addition, our model 
neglects corrections caused by quantum mechanical tun- 
neling. The classical avoided crossings are accompanied 
by quantum mechanical anticrossings, causing additional 
curvatures for intcrdot tunnel splittings comparable to 
the electrostatic interdot coupling energies^ 



IV. CHARGE AND TRANSPORT 
MEASUREMENTS 

In this section we discuss the measured ground state 
stability diagram of the TQD structure in Fig. [T] In the 
data presented we always use the left QPC (see Fig. [1]) 
as charge detector^ The differential conductance of the 
TQD is investigated within a range of the stability dia- 
gram that allows co-tunneling at strong enough tunnel 
couplings. 

The stability diagram in Fig. [4] displays the transcon- 
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Figure 4: (Color online) Measured charge stability diagram of 
the TQD device shown in Fig. [1] as a function of gate voltages 
V a and Vy. The color scale measures the transconductance of 
the left QPC in Fig.[T]as a function of the (modulated) plunger 
gate voltage Vy. Stable charge configurations are denoted 
by triples of numbers (Na, Nb, Nc). The graph is composed 
from several consecutive measurements explaining e. g. the 
horizontal line at V a ~ —0.22 V. 



ductancc dlqpc/dV^ of the left QPC as a function of 
gate voltages V a and V n with constant voltages applied 
to all other gates. Fig. [4] clearly demonstrates that a sin- 
gle nearby QPC is sensitive enough to monitor charging 
events in all three QDs. To measure transconductance 
using a lock-in amplifier, V 7 is modulated with an ampli- 
tude of AVy ~ 0.7 mV at a frequency of / = 33 Hz. In 
all transconductance measurements shown in this article 
V\ = — 300 //V is applied to ohmic contact I (compare 
Fig. [1]) in order to bias the QPC, if not stated otherwise. 
Measurements with a smaller bias voltage applied to the 
QPC assure that the linear response condition is fulfilled 
at V\ = —300 /xV. All other ohmic contacts besides con- 
tact I are grounded. The data in Fig.Q]fcatures lines with 
three different main slopes as expected for a TQD (com- 
pare with Fig. [3]). Almost horizontal lines of minimum 
transconductance are charging lines of QD A. Likewise, 
almost vertical charging lines belong to QD C. The slope 
of dV a /dVy ~ — 1 of the third kind of charging lines be- 
longs to QD B and shows, that the electrostatic coupling 
strengths between QD B and the two plunger gates a 
and 7 have similar values. The absence of all three kinds 
of charging lines in the lower left corner of the stability 
diagram suggests that the TQD is completely uncharged 
in this area of Fig. [4j 

However, the tunnel barriers of the QDs are larger for 
smaller gate voltages. In principle, charging lines can 
be invisible at very high tunnel barriers, if the charging 
process of a QD is slow compared to the time scales lim- 
iting the experiment. We ruled out the slow tunneling 
rate scenario by conducting careful test measurements, 
including different voltages applied to other gates than 
a and 7. We conclude therefore, that our TQD is really 
uncharged in the lower left corner of Fig. [4] 

Figure [5] expands a region of the stability diagram 
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Figure 5: (Color online) Expansion of a region of the TQD 
charge stability diagram in Fig.|4]for similar surface gate volt- 
ages. A black rectangle marks an area also marked by rect- 
angles in Fig. [7] 

in Fig. [4] for similar gate voltages applied. The charge 
reconfiguration lines (of positive transconductance and 
positive slope) are well resolved. The electrostatic in- 
tcrdot coupling (proportional to the length of charge re- 
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configuration lines) between the distant QDs A and C is 
small compared to those between neighboring QDs. In 
addition, the intcrdot coupling between QDs A and B 
is larger than that between QDs C and B. This is re- 
flected in the length of the charge reconfiguration lines 
between configurations (0,1,3) <-► (1,1,2) compared 
to (0,2,3) <-» (1,1,3). For few electrons charging the 
TQD, we find electrostatic interdot coupling energies of 
E AB ^ 680 /xeV, E BC ~ 150 fieV and E AC ~ 70A*eV. 
From the distances between charging lines we find charg- 
ing energies of the order E A ~ 1.1 meV, E B ~ 2.0 me V 
and Eq — 1.0 meV. The conversion of gate voltages to 
energies is done with the help of non-linear transport 
measurements as will be discussed at the end of this sec- 
tion (compare Fig.[7|). 

The model stability diagram in Fig. [3] was calculated 
for the energies and capacitances derived from the mea- 
sured stability diagrams. We find a good agreement of 
the main features, including the mean distances between 
charging lines and triple points and the average slopes. 
A complete quantitative agreement is not expected, be- 
cause of the limits of the model as discussed in section 
III. Moreover, the measured data reveal a spectrum of 
phenomena, not accounted for in the simple electrostatic 
model assuming a constant capacitance matrix. Some of 
these features are discussed below. 

The interdot coupling between QDs A and C increases 
as Vj is increased (compare the charge reconfiguration 
lines in Fig. \5\ between configurations (0, 1, 3) <-> (1, 1, 2) 
with (0,2,4) <-> (1,2,3) and (0,2,5) <-> (1,2,4)). This 
effect, that can also be observed when increasing V a , 
can be explained by considering two aspects. First, the 
quantum mechanical tunnel coupling between neighbor- 
ing QDs increases with increasing gate voltages, adding 
to the classical avoided crossing and causing a curvature 
of charging lines at the triple points. Second, the charge 
distribution and the position of the center of charges 
within all three QDs depend on the charge configuration. 

The distance between the almost horizontal charging 
lines of QD A varies strongly. A detailed analysis yields a 
charging energy of QD A that is larger for the third elec- 
tron than for the first or fourth electron (compare Fig.[2|. 
Such a strong effect implies a quite asymmetric confine- 
ment potential of QD A for the gate voltages applied. 
For instance, a situation where the first two electrons fit 
next to each other into an elongated QD, could explain 
the observation. 

Charging lines belonging to different QDs have differ- 
ent brightness, reflecting the amplitude of the transcon- 
ductance extrema. The brightness of a charging line is a 
linear function of its slope, diminishes at a larger width 
and is proportional to the electrostatic coupling between 
the QD and the detector QPC. The slope of a charging 
line determines the component of the line width parallel 
to the gate voltage Vy, which is modulated and relevant 
for the transconductance measurement. The quantum 
life time of an electron in a QD directly influences the 
widths of the charging lines. Generally, charging lines 



widen as gate voltages are increased and the involved tun- 
nel couplings grow. The charging line of QD B is bright- 
est, because QD B is well decoupled from the leads and 
exhibits the longest quantum life time. The widths and 
brightness of a charge reconfiguration line depends on the 
corresponding interdot coupling rather than a quantum 
life time. 

In the upper right corners of Figs. 2] and [5] charging 
lines turn into two parallel double lines, one with a large 
negative and one with a large positive transconductance. 
They are caused by current flowing through the TQD 
and the grounded contact III, branching off the current 
flowing from the biased contact I to the grounded con- 
tact II (compare Fig. [JJ). Note, that the ohmic contacts 
have resistances in the order of a R ~ 500 fl, resulting in 
a small potential drop across the TQD. A current max- 
imum at contact III causes a dip of the measured cur- 
rent at contact II. Hence, the transconductance dlu/dVy 
splits in a negative and a positive contribution, as ob- 
served. Comparison with Fig. [7] confirms the areas of 
enhanced transport. A finite current through the TQD 
along charging lines is caused by higher order tunneling 
processes. A detailed discussion of transport through the 
TQD follows at the end of this section. 

Within the triangular area marked with (0, 0, 1/0, 1, 1) 
in Fig. [51 the TQD fluctuates between the two charge 
configurations. We expect this bistability to be generic 
for serial systems of more than two QDs in the limit 
of small tunnel rates as will be discussed in a separate 
publication. 

Along the dark vertical line visible in Fig. O at V 1 ~ 
—0.272 V, the extrema of the transconductance appear 
to be more pronounced. This is caused by an internal 
switching in one of the measurement instruments and is 
not related to the TQD. 

Fig. [5] demonstrates the three dimensional nature of 
the TQD charge stability diagram. It shows a two- 
dimensional slice spanned by gate voltages Vp and V 7 
perpendicular to that in Fig. O which is spanned by V a 
and V 1 . The two stability diagrams in Figs. [5] and [5] to- 
gether allow to determine all capacitances necessary for 
model calculations with eq. [2] As before, charging lines 
with three different slopes can be identified, but two of 
the charging lines have similar slopes. They indicate a 
comparable ratio of couplings to both plunger gates /3 
and 7 for the two QDs A and B. In addition, as will be 
discussed in section V, the data in Fig.[B]contain features 
characteristic for the regime in which all three QDs are 
nearly energetically degenerate. 

Figure [JJ shows the differential conductance 
d^TQD / dVsD of the TQD as a function of the plunger 
gate voltages V a and V 7 for various bias voltages 
— 0.3mV < Vsd < 0.5 mV between the source and 
drain contacts (II and III in Fig. QJ of the TQD. The 
differential conductance is measured by means of lock-in 
technique with an ac modulation of AVsd = 20 fxV at a 
frequency of / = 33 Hz. For better comparison all gate 
voltages are identical for the transport measurements 
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Figure 6: (Color online) TQD charge stability diagram as in 
Fig. but as a function of the plunger gate voltages Vp and 
V-y . The graph is composed out of three independent measure- 
ments. The white region was not investigated. Slight changes 
of the internal potential between the measurements lead to 
slightly imperfect seams between parts of the graph (e. g. at 
V ~ -0.37 V; Vj ~ -0.28 V). Two dashed lines are a guide 
to the eye and follow very faint charge reconfiguration lines 
of positive transconductance. The bistable region observed in 
Fig. [S]is also visible. 

shown in Fig. [7] and the charge detection measurement 
displayed in Fig. [5] Note that, compared to Fig. [5j the 
area spanned by V a and V 7 is smaller for the transport 
measurements in Fig. [7J It corresponds to the upper 
right corner of the stability diagram in Fig. [5J where 
the tunnel couplings of the TQD arc largest. The 
logarithmic color scale for the differential conductance 
in Fig. [7J overemphasizes very small currents through the 
TQD. 

The differential conductance in linear response for 
Vsd — (Fig. [7)}) is entirely caused by higher order tun- 
neling processes as this area of the stability diagram does 
not contain quadruple points (compare discussion in sec- 
tion V). Alternatively, such higher order processes could 
be explained in a picture using molecular eigenstates^ 
of a tri-atom. Here, we restrict ourself to the picture of 
higher order tunneling processes in the basis of single-dot 
cigenstates. 

As the plunger gate voltages in Fig. [Jji are increased 
the tunnel barriers decrease and hence, the differential 
conductance increases. An exeption to this rule can be 
seen along the charging lines of slope dVa/dVj ~ — 1 
to the left of the black square in Fig. [7^ , along which 
the central QD B is in resonance with the 2DES in the 
leads. Here, the differential conductance is larger than 
along charging lines in the direct vicinity of the black 
square (for larger Vy). Near the black square only QD A 
or C can be resonant with the chemical potentials of the 
leads. Here, if away from triple points, current through 
the TQD is caused by third order tunneling processes 
through two non-resonant QDs in series. In contrast, 
along the charging lines of QD B an electron can occupy 
a resonant state in QD B between two sequential co- 
tunneling processes. Transport via such two successive 
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Figure 7: (Color online) The differential conductance 
d^TQD /dVsD measured through the TQD plotted with a log- 
arithmic color scale as a function of the plunger gate voltages 
V a and V~f. Voltages applied to other gates are as for Fig. [5] 
Black rectangles mark the same region as in Fig. [S] White 
color indicates a differential conductance exceeding the full 
range of the amplifier, small areas of the darkest blue color 
denote negative differential conductance. Applied bias volt- 
ages Vsd are indicated. 



second order processes is highly enhanced compared to 
transport involving one resonant first order and one third 
order tunneling process. 39 

For Vsd ^ the charging lines in Fig. [7] split into two 
parallel lines corresponding to two different resonances 
of one of the QDs with the source or drain potential in 
the two leads. For the same reason triple points turn 
into triangles. Similar as for a DQD^ the size of such 
a triangle or the distance between parallel double lines 
can be used to determine the conversion factors between 
gate voltages and the energy scales. These conversion 
factors are needed in order to calibrate the ground state 
stability diagrams and determine the charging energies 
and electrostatic interdot coupling energies of the QDs. 



V. VICINITY OF QUADRUPLE POINTS 

Because of the presence of three discrete charges in 
the triple-dot system, the charge-configuration diagrams 
are formally defined in a three-dimensional phase space. 
By analogy with the two-dimensional honeycomb dia- 
grams seen for double-dot structures, we term this three- 
dimensional charge configuration diagram a beehive dia- 
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Figure 8: (Color online) Three-dimensional charging dia- 
gram, or beehive diagram, showing stable configurations of 
the triple-dot system as a function of the three plunger-gate 
biases, normalised by their capacitive couplings to the closest 
dot. The front and bottom region of the structure corresponds 
to the configuration (0,0,0). On each of three visible end- 
faces resembles a two-dimensional honeycomb-like diagram, 
although more complicated diagrams can be seen at other 
slices, as shown in Fig. [5] 



gram, and a calculated example is shown in Fig. [5] In this 
case, we have deliberately chosen a capacitance matrix 
regime with minimal cross-coupling, so that the planes of 
the three visible end-faces show honeycomb-like charging 
diagrams. 

In the remainder of this article we focus on an area of 
the stability diagram where all three QDs of the TQD 
are close to being resonant with the chemical potential 
in the adjacent 2DES, that is where charging lines of all 
three QDs are close by. 

In the case of a DQD avoided crossings of any two 
charging lines always result in two triple points enclosed 
by four different areas of stable charge configurations, 
since two QDs are charged each by up to one additional 
electron (2 2 =4). In a TQD charging lines of three 
different slopes (belonging to the three QDs) exist. If 
two of them meet in a two-dimensional stability diagram 
they form triple points just as it is the case for a DQD 
(see Fig. [3]). In a three-dimensional stability diagram 
of a TQD, e. g. spanned by the plunger gates a, (3, and 
7, charging lines turn into planes and triple points turn 
into lines, (compare Fig. [5]). In a region where charging 
planes of all three QDs meet, each QD can be charged by 
one additional electron. This results in 2 3 = 8 possible 
charge configurations, surrounding four avoided crossings 
with eight triple lines. Such a three dimensional struc- 



ture contains four quadruple points, where two charging 
planes and two charge reconfiguration planes meet. Only 
at these quadruple points is transport by sequential res- 
onant tunneling of electrons through a serial TQD possi- 
ble. 

To more clearly examine the three-dimensional stabil- 
ity diagram in an experimentally accesible fashion, we 
study a series of parallel two-dimensional slices through 
the beehive diagram. We concentrate on regions of high 
degeneracy and use the terms appropriate for two di- 
mensions as triple points and charging lines. In Fig. [5] 
such measurements are plotted as a function of V a and 
V~t and compared with model calculations. The voltage 
ranges of V a and V 1 are identical for all subplots. The 
third plunger gate voltage Vp is increased in steps of 
2mV between -396 mV < Vg < -390 mV from (a) to 
(d) . The two bottom rows of Fig. [9] show conductance 
measurements, that will be discussed later. The middle 
row plots the transconductance of QPC charge detection 
measurements of comparable regions of the stability di- 
agram. The two upper rows feature identical numerical 
calculations according to our model, displayed with two 
different methods. 

The transconductance is measured with the left QPC 
with a bias voltage of Vu = —300 /iV applied to con- 
tact II (compare Fig. [T]). This bias voltage also causes 
current through the TQD at certain places of the sta- 
bility diagram. This is proven by the conductance mea- 
surements in linear response plotted in the two bottom 
rows of Fig. [5J Thus, the finite bias applied to contact 
II generates additional features in the transconductance 
measurements. These include extra lines, e. g. a line 
with slope —1 within the (2,2,3) region in Fig. O(biii), 
and gaps that interrupt lines, e. g. on the bottom left of 
Fig. O(diii). A small shift between the position of fea- 
tures in the transconductance data as compared to the 
conductance data can partly be explained by the applied 
biases, but could as well be caused by potential drifts 
during the time gap between these experiments. For a 
rough compensation the conductance measurements in 
the lowest two rows of Fig. [5] are horizontally shifted by 
AVg = 1 mV. Unfortunately, no charge detection mea- 
surements with the TQD left unbiased exist so far. Nev- 
ertheless, all features of the model calculations (first two 
rows of Fig. [5]) arc clearly seen in the measured transcon- 
ductance data (third row of Fig. [5]). 

Charge configurations, identified from a larger area 
stability diagram, arc depicted in the first row of the 
model stability diagrams of Fig. [91 where configurations 
X and Y equal (1, 2, 3) and (2, 1, 4), respectively. Placing 
the stability diagrams from Fig. [9] (a) to (d) above each 
other, with distances corresponding to Vg, results in a 
three-dimensional section of the stability diagram. Re- 
gions X and Y are similar to irregular pentangular-based 
pyramids inverted with respect to each other, (as can be 
seen in Fig. [5]). The tips of the pyramids arc oriented in 
approximately opposite direction from each other at two 
quadruple points. 



-0.1 f 



> 
-0.24 





,.Y ^^^^ 






(aii)| 


|(bii|i| 




piil 




m 

p 













- 0.06 



0.03 




10" 
10 



—6 h- " — " 



-0.2 V y (V) -0.1 E 



Figure 9: (Color online) Expansion of a region of the stability diagram of the TQD as a function of the plunger gate voltages 
V a and V-y. Here, charging lines of all three QDs are close by. Vp is increased in steps of 2mV from (a) to (d) or from (e) to (i). 
The upper two rows show identical results from model calculations in two different representations. Triples of numbers denote 
stable charge configurations, where X = 1, 2, 3 and Y = 2, 1, 4. The third row shows transconductance measurements as e. g. in 
Fig. [5] The two lowest rows display identical differential conductance measurements of the TQD with a linear (upper) and a 



logarithmic (lowest row) color scale (see main text for more explanations), 
all subplots and shown in the left corner. 



The voltage ranges of V a and V 1 are identical for 



For the approximately symmetric case shown in 
Fig-E^b), the central charging line with slope dV a /dV 1 ~ 
— 1 (belonging to QD B) shows a zig-zag behavior re- 
sulting in four triple points. In addition, the other two 
charging lines both contain a short segment parallel to 
the central charging line with dVa/dV^ ~ — 1. These fea- 
tures can be roughly explained as follows: Imagine the 
central charging line would be absent. Then we were 
left with one avoided crossing, where four lines end in 
two triple points. The central line, once added, repels 
the other four charging lines and four new avoided cross- 
ings occur, resulting in the observed geometry with eight 
triple points. 



labels in Fig.[5](bi)) 



In the following discussion we use a notation that sub- 
tracts the common charge state (1, 1,3) to be left with 
configurations of type (it, v, w) with u,v,w = 0, 1. In the 
most symmetric case the eight triple points then involve 
the following degenerate charge configurations (compare 



TP A 


(0,0,0)- 


+ (1,0,0)* 


+ (0,1,0) 


TP b 


(0,1,0)- 


-> (0,0,1)- 


- (0,0,0) 


TP C 


(1,0,1)- 


+ (1,0,0)- 


+ (0,1,0) 


TP D 


(0,1,1)- 


+ (0,1,0)- 


+ (0,0,1) 


TP e 


(1,1,0)- 


+ (1,0,1)- 


+ (1,0,0) 


TP F 


(1,0,1)- 


+ (o,i,i)- 


+ (0,1,0) 


TP G 


(1,1,1) - 


+ (1,1,0)- 


-(1,0,1) 


TP h 


(1,0,1)- 


+ (o,M) + 


-(1,1,1) 



At TPa either QD A or B can be resonantly occupied 
by an additional electron from the leads, but the occupa- 
tion of QD C is energetically forbidden. Hence, at TPa 
sequential tunneling of an electron through the TQD re- 
quires one co-tunneling process via an energetically for- 
bidden state in QD C. Similarly, transport at any other 
triple point requires one second order tunneling process. 
While TPa and TPb allow sequential co-tunneling of an 
electron, TPq and TPh allow sequential co-tunneling of 
a hole. Second order transport through the other four 
triple points involves two particles. 

Where two charge reconfiguration lines (blue) with 
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positive slopes and positive transconductance meet, two 
triple points combine to a quadruple point. As a func- 
tion of Vp this is possible for triple points TPa and TPb 
(close to the situation in Fig. [5](a)), triple points TPq 
and TP H (between Fig.[9](c) and (d)), triple points TP C 
and TPe (close to the situation in Fig. 0(c)), and triple 
points TP D and TP F (between Fig. 0(a) and (b)). The 
resulting quadruple points involve the following degener- 
ate charge configurations 

QP AB : (0,0,0)^(1,0,0)^(0,1,0)^(0,0,1) 

QP CE : (1,0,0)^(0,1,0)^(1,1,0)^(1,0,1) 

QP DF : (0,0,1)^(1,0,1)^(0,1,1)^(0,1,0) 

QP GH : (1,1,1)~ (1,1, 0)^(1,0,1)^(0, 1,1). 

At these four quadruple points resonant tunneling 
through the TQD is possible, e. g. at QPab an electron 
can sequentially tunnel from the left lead into QD A, QD 
B, QD C, and then escape to the right lead (or vice versa). 
Quadruple points Q AB and Qqh, respectively, allow se- 
quential tunneling of an electron versus a hole through 
the TQD, similar to triple points in a DQD. However, the 
nature of transport at quadruple points Q DF and Q CE 
can not be described by one electron or hole tunneling 
through the TQD, but involves two particles (electrons 
or holes). This extends the possibilities in a DQD, where 
only electron- or hole-like transport is possible. 

As a function of Vp the pairs of triple points TPc and 
TPe as well as TPd and TP F meet in the corresponding 
quadruple points, respectively, and then diverge again. 
During this process at the quadruple point one resonant 
charge configuration is exchanged between a pair of triple 
points, resulting in the modified triple points 

TP C : (1,0,0)^(0, 1,0)^(1,1,0) 
TP E : (0,1,0)^(1, 1,0)^(1,0,1) 

TP^ : (0,1,0)^(0,0,1)^(1,0,1) 
TP F : (0,0,1)^(1,0,1)^(0,1,1). 

Compared to the approximately asymmetric case in 
Fig.Etb) in Fig. Eta) TP D and TP F are replaced by TP' D 
and TP F and in Fig. (Sfd) TP C and TP E arc replaced by 
TP' C and TP E . 

Pairs of quadruple points as QPab Q^gh as 

well as QPq E *-* QPdf show electron-hole symmetry, 
respectively. The same is true for triple points, e. g. 
TPa «-> TPh- In addition, after subtraction of the com- 
mon charge state (1, 1, 3), triple points are pairwise point 
symmetric in respect to the central QD B regarding their 
charge occupation, e.g. TPa <-> TPb- 

The regions X and Y of stable charge configurations 
read X = (0,1,0) and Y = (1,0,1) after subtraction of 
the common charge state (1,1,3). Crossing the line of 
minimum transconductance, separating these two areas, 
from X towards Y involves adding a charge to QD A (or 
C). However, this is only possible via a QCA-process, 



where simultaneously one electron is pushed from the 
central QD B into QD C (or A). This is a combination of 
charging one QD and a charge reconfiguration between 
the other two QDs. Therefore, the slope of the QCA-line 
between regions X and Y is determined by the combi- 
nation of the two processes involved and differs from all 
other charging line slopes in the stability diagram. 

Let us now consider the reverse process, which involves 
crossing the QCA-line from Y towards X. During this sec- 
ond order tunneling process an electron leaves QD C (or 
A) and simultaneously pulls another electron from QD 
A (or C) into the central QD B. Interestingly, the com- 
bination of both processes (crossing the QCA-line for- 
and backwards) can result into transport of one electron 
through the TQD via two successive second order tun- 
neling processes, similar as along the charging line of QD 
B. 

Second order tunneling processes that preserve charge 
are usually called co-tunneling processes. The QCA pro- 
cesses described above are not charge preserving, but sec- 
ond order. Hence, we refer to these as QCA- co-tunneling 
processes. 

Note, that an equivalent situation to that shown in 
Fig. [5] occurs in the upper right quarter of the stability 
diagram in Fig. [6l but here, as a function of the plunger 
gate voltages Vp and . In Fig. [5] two charge transfer 
lines are retraced by dashed lines as a guide to the eyes. 
Because QDs A and B feature comparable electrostatic 
couplings to both plunger gates j3 and 7, some of the 
triple points are hardly seen in Fig. [5] (see also above 
discussion of Fig. [6|). 

The lowest two rows of Fig. [9] display the conduc- 
tance of the TQD plotted both with a linear (second 
lowest row) and with a logarithmic (bottom row) color 
scale. The conductance is measured in the linear re- 
sponse regime and for zero dc bias on all ohmic contacts. 
Comparison of the logarithmic conductance representa- 
tion with the model calculations in Fig. [9] shows, that 
along the charging lines belonging to the central QD B 
and at the QCA-line a small current flows through the 
TQD. Both kind of charging lines are distinguished, be- 
cause they allow transport through the TQD via two suc- 
cessive second order tunneling processes. As discussed 
above, along the other charging lines of QDs A and C 
only third order tunneling processes can cause transport. 
Accordingly, no current can be observed along charging 
lines belonging to QDs A and C. 

The linear representation of the conductance through 
the TQD (second lowest row in Fig. [9]) reveals distinct 
current maxima at quadruple points, and triple points 
near quadruple points. A detailed comparison with the 
model calculations in Fig. [9] suggests, that the conduc- 
tance maxima in Fig. (Uci) and [9]ji) arc very close to 
the quadruple points QPab an d QP GH , respectively. 
Fig. [9^6) and [UJgi) each show four bright maxima. The 
lower left one in Fig. [9{fi) corresponds to the quadruple 
point also seen in Fig.^ei). The upper left two maxima 
mark the triple points TPc and TP E in close vicinity of 
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QP CE , and the lower right maximum is close to QP^f- 
Figure Efgi) can be described accordingly. 

Strikingly, the QCA-line connecting TPc and TPp in 
Fig. [Ubi) is also visible as a line of minimal transcon- 
ductance in Fig. [UJfi) and [9jgi). On the other hand, 
the current flowing at the charging line of QD B is too 
small to be seen in the linear representation. Both trans- 
port channels involve two successive second order tunnel- 
ing processes. Still, the QCA-line near quadruple points 
shows a larger conductance than the charging lines of 
the central QD B. This suggests that QCA-co-tunneling 
processes, with two particles simultaneously moving, re- 
sult in a larger tunneling probability than regular second 
order co-tunneling processes, that can be explained in a 
one-particle picture. The origin of this phenomenon lies 
in the electrostatic interaction between two electrons tun- 
neling simultaneously and will be discussed in a separate 
publication. 

Finally, we would like to note that spin blockade of 
transport in linear response through a TQD is expected 
for certain quadruple point configurations. It was not 
observed in the regime treated in Fig. [5J This can in 
part be explained by the special configuration, where up 
to three electrons are added to the charge state (1, 1,3), 
but at most one electron to each QD. In the configuration 
(1, 1,3), each QD already has a spin 1/2. After an extra 
electron charge has been added to one of the QDs, this 
QD has zero spin. This implies that this QD can now 
provide an electron with arbitrary spin (up or down) to 
tunnel to an adjacent QD. Hence, a full spin blockade 
is not expected for the region of the stability diagram 
discussed in Fig. [9] 

VI. CONCLUSION 

In summary, we have realized a lithographically de- 
fined serial triple quantum dot that can be tuned to con- 
tain any number of electrons between zero and about 



ten in various configurations. Quantum point contacts 
as integrated charge sensors allow to determine the ex- 
act number of electrons charging each of the quantum 
dots. We have studied the ground-state stability diagram 
of the triple quantum dot in close vicinity to quadru- 
ple points where four different charge configurations are 
energetically degenerate. In this regime, quantum cel- 
lular automata processes are observed among other fea- 
tures, adding to the physics that can be found in dou- 
ble quantum dots. A simple electrostatic model, that 
can easily be scaled to larger structures containing more 
than three quantum dots, is compared with our data. 
A detailed discussion of the conductance near quadruple 
points reveals several kinds of tunneling processes. Quan- 
tum cellular automata co-tunncling processes lead to an 
enhanced conductance at only two-fold degeneracy of the 
triple quantum dot. The excellent control of charge con- 
figurations and tunnel couplings achieved in this triple 
quantum dot now opens the possibility to study coher- 
ent dynamics, e. g. of charge and possible spin transfer in 
such a complex quantum system. 
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